Production of magnetic energy by macroscopic turbulence in GRB afterglows 
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ABSTRACT 

Afterglows of gamma-ray bursts are believed to require magnetic fields much stronger 
than that of the compressed pre-shock medium. As an alternative to microscopic plasma 
instabilities, we propose amplification of the field by macroscopic turbulence excited by 
the interaction of the shock with a clumpy pre-shock medium, for example a stellar wind. 
Using a recently developed formalism for localized perturbations to an ultra-relativistic 
shock, we derive constraints on the lengthscale, amplitude, and volume filling factor 
of density clumps required to produce a given magnetic energy fraction within the 
expansion time of the shock, assuming that the energy in the field achieves equipartion 
with the turbulence. Stronger and smaller-scale inhomogeneities are required for larger 
shock Lorentz factors. Hence it is likely that the magnetic energy fraction evolves as the 
shock slows. This could be detected by monitoring the synchrotron cooling frequency 
if the radial density profile ahead of the shock, smoothed over clumps, is known. 



1. Introduction 

Since their discovery, gamma-r ay burst afterglows hav e been attributed to synchrotron radi- 
ation from the forward sho c k wave (IMeszaros &: ReesI 119971 ). although it has been recently argued 
(jUhm &: Beloborodovl 120071 : iGenet et al.l 120071 ) that observations might support a model in which 
the forward shock is invisible and the afterglow is emitted by a long-lived reverse shock in the burst 
ejecta. Assuming anyway a forward-shock origin for the afterglow emission, it is difficult to account 
for the magnetic energy density behind the forward shock by simple compression of the pre-shock 
field. Interstellar magnetic energy densities are typically comparable to thermal pressures and are 
therefore a fraction £3,0 = Pmag,o/po = — of the total internal energy density when rest 
mass is included. It is possible that the pre-shock medium is the stellar wind of the burst progeni- 
tor; while the magnetic energy fraction in winds is less well known, it is unlikely to be much larger 
than this. Simply compressing the medium would produce approximately the same ratio Pma.g/ P 
behind the shock. Instead, phenomenologica l models of after g low light curves typically require 
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In fact, gamma-ray burst (hereafter GRB) afterglows present the most compeUing case among 
astrophysical coUisionless shocks for prompt creation of magnetic energy. The synchrotron emission 
from supernova remnants is generally consistent with compression of the int erstellar field, alth ough 
some modest additional amplification may be required in particular cases (|Volk et al.ll2005l ). We 
focus here on GRB afterglows rather than GRB internal shocks, which also emit by some com- 
bination of synchrotron and synchrotron self-Compton, because it ma y be that the burst ejecta 
are magnetically dominated (ICoburn &: Bogga l2003l : IZhang et alJ l2003l ) , and so the internal shock 
emission is driven by the progenitor magnetic field. 

The leading hypothesis for field amplification in GRB afterglows is the relativistic Weibel 
instability, which extracts free energy from the anisotropy of the particles' velocity distribution 
function, producing filamentary currents aligned wi th the shock normal; the se currents are respon- 
sible for the creation of transverse magnetic fields (iMedvedev &: Loeblll999l ). This process is able 
to violate MHD flux-freezing because it occurs on a microscopic scale — the relativistic electron or 
ion skin depth — where the inertia of individual charged particles is significant. While the Weibel 
instability provides a plausible mechanism to isotropize the particle velocities, it is unclear whether 
the small-scale fields that it produces can survive mutual annihilation long enough to explain the 



observed synchro tron afterglow emission. Several groups (ISilva et al 



2003 



Frederiksen et al 



2004 



Spitkovskvl I2OO5I ) have attempted to simulate the long-ter m nonlinear ou tcome of the instability, 
but a consensus on this question has not been achieved yet (|Waxmanll2006l . and references therein). 
One might have thought that if this instability were the source of post-shock fields, then ea should 
have a universal value for highly relativistic, highly coUisionless shocks. Yet, while is mod- 
eled by a constant for indiv i dual GRB afterglows, it seems to vary from one afterglow to another 
(jPanaitescu fc Kumail[2002 : IVost etaP 1200.4 IPanaitescJ[2005l ). 



In this paper, we explore a traditional magnetohydrodynamic explanation for magnetic field 
growth: turbulence. It is well known in non-relatiy istic fluid dynamics that oblique shocks produce 
or alter the vorticity of a fluid (jlshizuka et al.lll964l ). In this paper we will show that the same is true 
for an ultra-relativistic shock passing over density inhom ogeneities in the pre-shock circ umburst 
medium. The formalism described in a previous paper (| Goodman &: MacFadyenI 120071 ) has let 
us deflne the vorticity created in an ultra-relativistic fluid in which the energy-momentum tensor 
can be approximated by that of an ideal fluid with pressure equal to one third of the proper 
energy density (P = p/3). In the same work, we have introduced a remarkably simple but accurate 
general approximation for the local modulation of the shock Lorentz factor (F) by pre-shock density 
inhomogeneities; within this approximation, it is not necessary to follow the details of the flow far 
downstream in order to predict the evolution of the shock, provided that F ^ 1, that the pre-shock 
pressure is negligible and that the post-shock pressure satisfles P = p/ 3. This approximation. 



which is modeled on non-relativistic results described bv IWhithamI (jl974l ). reproduces exactly the 
self-similar evolution of F for a shock advancing into a cold pre-shock medium with a power-law 
density proflle in planar, cylindrical, or spherical symmetry ()Sarill2006l ). More importantly for the 
present purpose, it allows us to estimate the post-shock vorticity resulting from a prescribed pre- 
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shock density that varies along as well as perpendicular to the shock normal. Given the vorticity, 
we divide the post-shock velocities, which are marginally non-relativistic in the average post-shock 
rest frame, into vortical and non-vortical parts. We presume that the energy density in vortical 
motions is a measure of the magnetic energy density that will eventually result after the eddies 
wind up the field to the point where its backreaction on the turbulence becomes important. These 
methods are described in §11. 

In §111, we briefly review the present state of knowledge concerning the density inhomogeneities 
that may exist ahead of the shock. Both the amplitude and the lengthscale of the inhomogeneities 
are important. The former controls the amount of vortical energy — and then of magnetic energy — 
that is produced, while the latter determines the eddy-turnover time of the turbulence, which — 
when multiplied by the number of eddy rotations necessary to amplify the field up to the observed 
€b — must be less than the shock deceleration time, so that the field can be significantly amplified 
before adiabatic expansion reduces the particle energies available to be radiated. The uncertainties 
are large because one doesn't know whether the pre-shock medium is more like a stellar wind or like 
some component of the Galactic interstellar medium, and because the lengthscales of interest are too 
small (< 10^^ cm) to be directly resolved even in the interstellar medium. Inhomogen eities on some- 



what larger scales have b e en invoked to explain undulations in aft erglow light curves (jWang &: Loeb 



2000 



Lazzati et al.ll2002l : ISchaefer et al.ll2003l : iNakar et al.l 120031 '). 



In §IV, we use the formalism of §11 to characterize the density contrasts and lengthscales that 
pre-shock clumps should have in order to amplify the magnetic field up to the observed value, in the 
light of the circumburst picture outlined in §111. We find that, for smaller shock Lorentz factors, the 
constraints on clump sizes and overdensities become less stringent; as a consequence, the magnetic 
energy fraction produced by pre-shock clumps via macroscopic turbulence is expected to evolve as 
the shock slows down. In §V, we comment on the plausibility of our proposed mechanism to explain 
the magnetic field amplification in GRB afterglows and we discuss how the results obtained in §IV 
could be tested by inferring the time dependence of from the time evolution of the observed 
synchrotron cooling frequency as the shock ages. 



2. Geometrical Shock Dynamics 



The evolution of a shock advancing into an inhomogeneous medium depends, in principle, upon 
the details of the downstream flow behind the shock and of the "piston" that drives it. Geometrical 
Shock Dynamics (hereafter GSD) is an approximation fo r this evolution in which only the conditions 
at the shock appear explicitl y. Originally formulated by I WhithamI (|1974| ) for non-relativistic fluids, 
GSD has been extended bv I Goodman &: MacFadvenI (|2007l . hereafter Paper I) to strong (F ^ 1) 
ultra-relativistic shocks advancing into an ideal fluid whose pressure is negligible ahead of the shock, 
but one third of its proper energy density behind the shock (P = p/3). 



The fundamental approximation of GSD is to evaluate the forward-going Riemann character- 
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istics in the post-shock flow as if that flow were (i) isentropic and (ii) homogeneous far downstream 
with properties determined by the mean shock speed and the mean pre-shock density. Actually, 
pre-shock density inhomogeneities lead to post-shock entropy variations, so assumption (i) is wrong 
in principle, but it turns out to be a useful fiction. In this respect, an ultra-relativistic flow has 
the advantage that since pressure depends only on energy density and not on any other thermody- 
namic variable (such as the proper number density of baryons, A^), the actual entropy is irrelevant 
to the Riemann characteristics, which therefore enjoy exact Riemann invariants. So assumption 
(i) is well justified except insofar as it may be compromised by secondary shocks created by the 
inhomogeneities themselves behind the main shock. Assumption (ii) is reasonable when pre-shock 
density inhomogeneities are small in lengthscale, so that they may be expected to average out far 
downstream. 

With assumption (i), the Riemann invariant on the forward characteristics has the same value 
just behind the shock as it does far downstream, and therefore, with assumption (ii), the same 
value that it would have in the mean flow. Together with the jump conditions across the shock, 
this provides a relation between the local Lorentz factor of the shock, F, and the proper pre-shock 
energy density />o ~ mN^c^, in which Nq is the proper number density of nucleons ahead of the shock 
and m is the rest mass per nucleon. In the essentially one-dimensional case that pQ varies along the 
shock normal but not perpendicular to it, the ultra-relativistic GSD relation for the response of the 
shock Lorentz factor F to localized and transitory variations in the pre-shock density po becomes 
(Paper I) 

F = f ( ^ ) where X = V^-\k 0.232 (1) 

\PoJ 2 

and the corresponding change in the post-shock pressure is 

/ \ 1-2A / X 0.536 

where the overbars indicate mean values. As in non-relativistic GSD, these relations can be extended 
to multidimensional flows in which pQ varies laterally as well as longitudinally (with respect to the 
shock normal), causing convergence or divergence of the shock normals. Relations ([T]) and ([2]) are 
then modified by factors involving the ratio of the local shock area to its mean value (Paper I). It 
is shown, however, that these corrections are of higher order in F^^ unless the density contrasts 
are ~ 0(F). For the conditions contemplated in this paper, eqs. ([T]) and ([2|) will be adequate even 
in two or three dimensions. 

As in the original non-relativistic theory, rigorous error estimates for ultra-relativistic GSD 
are difficult. Informally, the following conditions are probably necessary for the approximation to 
be useful. First, the pre-shock medium should be cold, meaning that pre-shock pressure satisfies 
Pq ^ Po ^-nd that internal velocities are <C c; this is very likely true of the external forward shocks of 
GRBs. Second, the lengthscales of the pre-shock inhomogeneities should be small compared to the 
shock radius, so that the shock responds to local perturbations before conditions far downstream 
have time to react. Third, since F~^ is used as a small parameter, the pre-shock density fluctuations 
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should not be so large as to cause the shock to become sub-relativistic, i.e. one requires po/po ^ 
fi/^^ Finally, transitions in density should not be so abrupt as to cause strong reverse shocks, 
which would alter the forward shock dynamics. Paper I describes tests of ultra-relativistic GSD 
by comparison with exact self-similar solutions (some of which it reproduces exactly) and with 
numerical simulations. The latter indicate, for example, that eqs. ([T|) and ^ are in error by only 
a few percent for a f = 10^ shock encountering overdensities as large as po/po < 30. 



2.1. Relativistic vorticity 

This subsection is independent of GSD. We review the meaning of enthalpy current and vortic- 
ity in ideal relativistic fluids, especially those with the ultra-relativistic equation of state P = p/3. 

When the shock passes over a local density excess — considered, for simplicity, in isolation 
from other inhomogeneities — the resulting post-shock velocities are of two kinds. First, since the 
shocked clump is overpressured compared to its post-shock surroundings [eq. ([2])], it will expand 
and drive an outgoing pressure wave. If the density contrast of the clump is small, then the wave 
is essentially a linear disturbance from the start and travels at the sound speed, c/\/3, in the 
rest frame of the mean post-shock flow; waves launched by large overdensities will be somewhat 
faster and may steepen into secondary shocks, but whatever its strength, the pressure wave rapidly 
departs its source. Overlapping pressure waves launched by many distant clumps may contribute 
significant local velocity perturbations, but, because of their oscillatory nature, intuition suggests 
that these velocities will not secularly amplify the magnetic field (except insofar as secondary shocks 
may contribute to vorticity — see §IV). It would be interesting to test this expectation in numerical 
simulations. 

Unless the density excess is constant along the shock front, the post-shock velocity field will 
also contain a vortical component, whose strength is estimated below for an initially spherical 
overdensity with a gaussian radial profile. As shown in Paper I, the equation of state P = p/3 
allows some freedom in how one defines relativistic vorticity, but to be useful in constraining the 
evolution of the flow, the vorticity should be associated with a conservation law such as Kelvin's 
Circulation Theorem, 

^^j>H^dx^ = Q, (3) 

c 

where C is a closed contour comoving with the fluid four- velocity U^, and = hU^ for an 
appropriate thermodynamic function h (see below). Equation ([3]) is equivalent to 

-^-Vx(t)Xcj) = 0, uj = cVxH, (4) 

where = cW/U^ is the three-velocity of the fluid, and H is the spatial part of H'^. Although 
formally identical to the non-relativistic vorticity equation and written with three- vectors, eq. ^ 
is actually relativistically covariant. 
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Because vorticity and circulation travel with the local flow velocity and are nonoscillatory 
in the local fluid frame, they have the potential to twist up and secularly amplify any magnetic 
fleld frozen into the flow. A fundamental assumption of this paper is that the magnetic energy will 
eventually reach equipartition with the kinetic energy invested in the vortical part of the flow. This 
assumption would also be well worth testing numerically. 

For a general equation of state P = P{p, N), where is the proper number density of conserved 
particles (e.g., baryons), it is conventional to take the quantity h in the relation H'^ = hU'^ to be the 
enthalpy per particle: h = {p + P)/N . One assumes an ideal fluid, so that the energy-momentum 
tensor is 

T^''' = {p + P)U>'U'' + Pg^"" , (5) 
{g^'^ ~^ diag(— 1, 1, 1, 1) in Minkowski coordinates), and the equations of motion are 

T^% = and (iVC/^);^ = 0. (6) 

The First Law of Thermodynamics dh = TdS-\-N^^dP implies h^^ = P^^/N \i the fluid is isentropic; 
the flrst equation of motion in eq. ^ can then be recast as 

U''H% = -h^, (7) 



the "curl" of which implies eqs. ([3]) and ([!]) (e.g.. lEshraghill2003l . and Paper I). These relations do 



not hold across shocks, of course, since there is an entropy jump. But if the pre-shock medium 
is inhomogeneous, then even if the flow behind the shock is smooth, it will not be isentropic in 
general; that is, P/N^/^ (for an ultra-relativistic equation of state) will not be uniform. Therefore, 
the conventional choice of h does not lead to a conserved circulation under the circumstances 
contemplated in this paper. Fortunately, as pointed out in Paper I, if one deflnes using h oc 
P^/^ instead of the true enthalpy, then eq. ([7]) always holds in smooth parts of the flow. This 
is a consequence of the equation of state P = p/3, which is "barytropic" if not isentropic. It is 
convenient to choose the constant of proportionality so that reduces to the fluid four-velocity 
when pressure is uniform. Therefore, we replace the conventional enthalpy current with 

Rf" = {P/pf^^Uf' . (8) 

With this choice, circulation is conserved [eqs. ([3]) & ^] everywhere except across shocks. 



2.2. Vorticity production by shocks 

The goal of this subsection is to use the one-dimensional GSD approximation to derive eqs. (jisp - 
()16p . which relate the post-shock vorticity to the pre-shock fractional overdensity 5 = (po/po) — 1- 
Figure [1] illustrates four stages in the interaction of an ultra-relativistic shock with a density clump. 

We begin by recalling some basic consequences of the shock jump conditions that will be 
needed below. Let [Q] denote the discontinuity in a fluid property Q across the shock front. In 



-7- 



the instantaneous local rest frame of the shock, where the outward unit normal to the shock front 
is n, the jump conditions are [T'^-^nj] = 0. Using eq. ([5]) for T^'^ (ideal fluid) and assuming that 
Pq Po ahead of the shock and P = p/3 behind it, one finds that the post-shock three-velocity of 
the fiuid is • n = — c/3 in the shock frame (hence subsonic, since the sound speed is c/Vs). The 
post-shock energy density is p = 2T^pq, where F ^ 1 is the local Lorentz factor of the shock in the 
rest frame of the pre-shock fluid. Similarly, conservation of particles implies [NWrii] = 0, whence 
N = 2\/^TNo. The quantities p, P, and N will always denote proper values, meaning that they 
are defined in the local fiuid rest frame and are therefore Lorentz invariants by fiat, with subscript 
"0" denoting a pre-shock value rather than a spacetime index. 

To facilitate Lorentz boosts between the pre-shock and post-shock or shock frames, it is often 
convenient to use the rapidity parameter tanh~^(f/c), where v is the three- velocity in the pre- 
shock frame. The relativistic addition of colinear three-velocities is equivalent to addition of the 
corresponding rapidity parameters. Thus, for example, using <I> = cosh~^ T ~ ln(2r) for the rapidity 
parameter of the shock and (p for the rapidity of the post-shock fiuid, it follows from the above that 
tanh((/) — <!') = —1/3, whence = $ — ln\/2. The Lorentz factor of the post-shock fiuid relative to 
the pre-shock frame is cosh(/) = cosh(^>-ln^/2) = r/V2 + 0{r-^). 

Pre-shock clumps will typically have comparable longitudinal and lateral dimensions (meaning: 
along and perpendicular to the mean direction of shock propagation) in their own rest frame. In the 
shock and post-shock frames, the clumps will be longitudinally contracted by factors ~ <C 1. 
During the transit of the shock over a clump, and even during the subsequent expansion of the 
shocked clump as it comes to pressure equilibrium with the surrounding post-shock fiuid, there 
will not be enough time for signals (sound waves) to communicate laterally from one end of the 
clump to the other. Therefore, the interaction of the shock with the clump can be calculated in a 
one-dimensional approximation, in which the area of the shock is constant and the pre-shock mass 
density (pq/c^) is stratified on planes parallel to the shock front. 

Nevertheless, lateral density gradients do produce post-shock vorticity, u, even in our one- 
dimensional approximation. First of all, the longitudinal component -ff y of the non-conventional 
enthalpy current defined in eq. ([8]) varies with lateral position behind the shock. Secondly, since the 
shock itself is delayed differently at different lateral positions, the shock normal develops a lateral 
component, which leads to a small lateral current, H±. Although the magnitude of H± is 0(F^^) 
compared to i^y for a given clump amplitude, its longitudinal derivative makes a contribution to 
uj that is comparable to the lateral derivative of i^y, as a consequence of Lorentz contraction. 

Let 5{r) be the density contrast of the clump: po{r) = pq[1 + S{r)], where po is the mean pre- 
shock value. The shock propagates along z on average; let us choose a cartesian coordinate system 
X, y, zq in the pre-shock rest frame, where the subscript "0" is used to distinguish between the post- 
shock (z) and pre-shock (zq) longitudinal coordinates. In our one-dimensional GSD approximation, 
it follows from eq. ([T|) that 

F(r) = f[l + 5(r)]-\ (9) 
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if f ^ 1 is the value for a smooth pre-shock medium. As discussed at the beginning of this section, 
it win be assumed that |(5|^ <C f so that F ^ 1 at ah times; since A ~ 0.232 is fairly small, rather 
hefty density contrasts can be compatible with this condition. 

First of all, let us focus on the contribution of H to the post-shock vorticity. Let T(r) be the 
time in the pre-shock frame at which the shock reaches pre-shock position r: 

cr(x, y, zo) ~ + j [l + <5(x, y, z'^)f^ dz^ = zo + ^f-2/(x, y, zq). (10) 

where we have assumed that the typical size of a density clump is small compared to the char- 
acteristic lengthscale for variations in F. The shock surface at time t is determined implicitly by 
T(r) = t, so that the shock normal is n = \W qt\~^W qt (the subscript reminds that derivatives 
are taken ahead of the shock). To leading order in f the lateral components of the normal are 
therefore 

20 

~ Af-^v^/ = Af-2 + 5'f^-\V^5') d^. (11) 

Here 5' is shorthand for 6{x,y,ZQ) and V_l = {d/dx,d/dy,0) is the same in the pre-shock and 
post-shock reference frames. The lateral part of the post-shock fluid 4- velocity is U± « nj_U^^ ~ 
n±T/^/2. With use of eqs. ([5]) and Q, the post-shock lateral enthalpy current becomes 

« 2'^/'^Xr'\l + (5)(i-^^)/^Vx/. (12) 

Of course the lateral components of the enthalpy current take the same values in the mean post- 
shock and pre-shock frames since these frames differ by a longitudinal boost. However, in order to 
compute the vorticity in the mean post-shock fluid frame, we should remember that the longitudinal 
derivative in this reference system Vy is related to the corresponding derivative Vo,|| in the pre- 
shock frame by Vy = cosh^ Vgji because of Lorentz contraction. Thus the contribution of H± to 
the post-shock vorticity is 

V|| X /fx ~ ^A(l + (5)(2^-=^)/^(n X V<5) 

+ ^^^^(l + <5)-3(i+2A)/4|i(^x VI) (13) 

o uZq 

where n = is the mean shock normal and clearly the differential operator n X V is the same in 
the pre-shock and post-shock mean rest frames. 

A comparable contribution to the post-shock vorticity comes from The longitudinal 

component of the post-shock 4- velocity is C/|| = sinh((/> — cp) when measured in the mean post- 
shock frame. Since <p — <j) = ^ — ^ = —A ln(l -|- 5), 

1 r. 



l + 5f-(l + Sr^ (14) 



in the mean post-shock frame. This tends to be quite subluminal: for unit overdensity {S = 1), for 
example, Un ~ — 0.1615 6;^, and the corresponding fluid Lorentz factor is ~ 1.013. The longitudinal 
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post-shock enthalpy current is = (P/P)^/^J7|| = Taking the lateral derivative 

of this and combining with eq. (jl3p and the definition of / in eq. (jlOp . we can write the total post- 
shock vorticity as 

+ c ~ (1 + ^)-3(i+2A)/4|i 5'f^'\n X V5') dz'^ , (15) 



dzQ 

where the coordinate system is in the pre-shock rest frame and the derivatives are taken ahead of 
the shock. We have introduced the function 



^(2A+3)/4 



6A -|- 1 ^, X 6A — 1 ^, X 
^— (1 + 6)^ + ^^(1 + S)-^ 



(16) 



For \5\ <C 1, f{6) PS 3A/2 0.348. It is interesting that in the total post-shock vorticity there are 
no surviving factors of f. 



2.3. Vortical energy 

We have not been able to formulate a rigorous relativistically covariant way of dividing the 
energy of the post-shock flow into vortical and non- vortical parts. In a non-relativistic flow, how- 
ever, this would be straightforward. One would divide the three-velocity field into potential and 
solenoidal parts, 

t, = VV- + V X A, (17) 

and then define the vortical energy by 

^vort = ^ y I V X A|2 d^X , (18) 

where pm is the non-relativistic mass density. To make the decomposition (jl7p unique, some mild 
additional restrictions are necessary: for example, that the region of interest is simply connected, 
and that ip and A have some specified behavior on the boundary or at infinity. Then one can 
impose V • A = and solve for A from 

V'^A = -uj, (19) 

with a; = V X I?, using an appropriate Green's function to invert V^. The nonvortical part follows 
similarly from V'^tp = V • i?. 

It isn't clear how to proceed in the relativistic case because the coordinate energy density T^^ 
isn't simply quadratic in v, in general. However, for the applications we have in mind, the vortical 
motions are plausibly subsonic in the mean post-shock frame (as long as 5 is not too large), and 
therefore only mildly relativistic, {v/c)"^ < 1/3. Therefore, we will use the decomposition (jl7p . 
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except for three changes: (i) cH replaces v, and therefore cj — > cV X ii" m eq. p^ : (ii) pm, the 
non-relativistic mass density, is replaced by p/c^, where p is the proper internal energy density; (iii) 
the factor 1/2 in eq. p!7|) is replaced by 4/3. The reason for (iii) is that if v/c <C 1 and P = p/3, 
then TOO p + (4/3) ^ (y/cf, and we assume that -Bvort derives from the second term0 

A pre-shock density clump will experience, after the transit across the shock front, a contraction 
of its longitudinal size by a factor q~^, where a = N/Nq ~ 2\/2r ^ 1. What was an approximately 
spherical clump becomes a pancake, with all of its associated vorticity in or near the interior. We 
adopt post-shock cylindrical coordinates R, 9, z with the z axis parallel to the mean direction of 
propagation of the shock. Taking the origin at the center of the clump, assuming the clump to be 
axisymmetric, and recalling that the vorticity computed in eq. (jl5p has been expressed in pre-shock 
coordinates, we may write u> = u){R, zq) sq = lj{R, az) eq a{R)5Y,{z)eQ, where 5j^{z) is the Dirac 
delta function, and 

+00 +00 
a{R)= j Lo{R,az)dz = ^ j u;{R,zo)dzo (20) 
—00 —00 
will be referred to as the "projected vorticity". 

Just after the shock passage, a clump will be at higher pressure than its surroundings, as already 
discussed, and will therefore expand until it reaches approximate pressure equilibrium. Since the 
sound-crossing time in the longitudinal direction will be much less than that in lateral directions, 
equilibrium will be reached with only a small fractional change in the lateral size. The fractional 
change in the longitudinal size will be of order unity, but for F S> 1 and clump overdensities 
6 <^ f'7ad/(i-2A) ^^^^ = 2 is the adiabatic index appropriate for the one-dimensional expansion of 
an ultra-relativistic fluid) the clump will remain highly flattened even after equilibration. Since the 
vortex lines lie parallel to the clump midplane (i.e., in the direction eg) and are "frozen" into the 
clump (as a consequence of Kelvin's Circulation Theorem), the expansion will have little effect on 
the projected vorticity in eq. (j20p . We may therefore solve eq. (jl9p in the approximation that the 
vorticity is confined to a thin sheet with the projected vorticity computed in eq. (j20p . The problem 
is mathematically equivalent to finding the vector potential due to an axisymmetric current sheet 
carrying a toroidal current: 

+00 

A{R, z) = A{R, z)eg = ^ j a{k) Ji{kR) e-'^^'Uk eg (21) 



^Why 4/3 rather than 1/2, as one might expect in the non-relativistic limit? The answer is that the rest- 
mass density of non-relativistic fluid mechanics is actually pm ~ jNm, where m is the mass per particle and 
7 = whereas the proper energy density is p = Nm{(? + u), where mu is the internal energy per particle. 
Thus, for a general equation of state, r°° = -f (7^ - 1)P = "/pmic^ + u) + (7^ - 1)P. For (v/c)^ < 1, this 
reduces to T"" ~ PmC^ + pm{v'^ /'2 + u) + [{v / cf' {pmu/2 + P)]. The middle term is now recognizable as the non- 
relativistic kinetic-plus-internal energy density. The term in square brackets would normally be negligible for a cold 
fluid because u ~ P/pm ^ c?\ however, for an ultra-relativistic ideal fluid u ~ 'iP'y/pm ^ (? and the result is 
r"" p^c? + p^v^ jl + 3P + 4P(?;/c)^, of which the first two terms are now negligible. 
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where Ji{kR) is the Bessel function of order 1 and a{k) is the Hankel transform of the projected 
vorticity a{R): 

+00 

d{k)= j a{R)Ji{kR)RdR . (22) 


As soon as the clump is at the same pressure as its surroundings, if the vortical motions are sub- 
relativistic we can use the decomposition (fTTl) with t?vort — > c vort = V X A. Moreover, pressure 
equilibrium between the clump and the average post-shock medium implies El that the clump proper 
energy density p is approximately equal (neglecting terms of order 0{v'^ /(?)) to the mean post-shock 
value p. Then the vortical energy of a single clump becomes 



E. 



vort 



In accordance with the discussion following (jl9p . we have replaced the factor 1/2 in (jl8|) with 4/3. 
Also, we have used integration by parts to replace the integral over all space with an integral over 
the source of vorticity only. This is particularly convenient when the source is represented by a 
vortex sheet, because we may express the energy as an integral in k space: 

+00 +00 
87r_ f A{R,0)a{R) Att _ f a^{k) 

EyoTt = —p / ^ RdR=—p / —^dk. (24) 



We shall assume that the vortical energies of different clumps can simply be added. This 
is justified if the clumps are well separated compared to their larger [i.e., lateral) dimensions. 
Then if the number density of clumps in the pre-shock frame is A^^o and all the clumps have the 
same axisymmetric density profile, the vortical energy density in the average post-shock frame is 
Pvort = ct ]^cfiEvorti and the vortical energy fraction becomes 

+00 

_ pvort 47r f a'^{k) 

Cvort = = ^aA^c,o / — 5— dfc. (25) 
p 3 J 


Recall from eq. ([20]) that the projected vorticity is proportional to = {2^/2t)^^ . Since this 
factor is squared in computing the vortical energy, it follows from eq. (I25p that Evort oc F^^ for a 
fixed pre-shock density field. This scaling is perhaps the most important conclusion of our analysis 
up to this point. 



^In the mean post-shock rest frame, if T'"' is the energy-momentum tensor of the fluid in turbulent motion 
and T''" for the mean fluid (whose four- velocity has spatial components i7i = in this reference frame), pressure 
equilibrium requires T'-' = T'-'; under the assumption of spatial statistical isotropy for the turbulent motions, this 
impUes T*'' = 3P, which yields, for an ultra-relativistic fluid (T'"'' = T^*" = 0), the equality = p; then, since 
T^" ~ p -f- (4/3) p (u/c)^, if the vortical motions are sub-relativistic we obtain p ~ p. 
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The circumburst medium 



Observations support the idea that long-duration GRBs are associated with the deaths of 
mass ive Wolf-Rayet (WR) stars, presumably arising from their core-collapse (jWooslev Sz Bloom 
20061 . and references therein). Then the circumburst environment is determined by the star's mass- 
loss history. At the onset of the WR phase, the WR stellar wind is expected to expand with 
a typical velocity fwa ~ 2000 kms~^ inside the pre-existing slower wind emitted during the red 
supergiant (RSG) phase, whose characteristic speed is Ursg ~ 20 kms~^. The win ds from massive 
RSGs are characterized by a mass-loss rate Mrsg between 10~^ and 10"'^ MQyr~^ ( Chevalier et al. 
2006 ) ■ while the mass-loss rates Mwr of WR stars are between 10^^ and 10"'^ MQyr""^ ( Crowthei 
20061 ^. Several solar masses are shed by the star during these evolutionary phases. The mass 
equivalent to the energy of a GRB, on the other hand, is only ~ 0.06£'iso,53 Mq. Therefore, 
the GRB forward shock is expected to become non-relativistic long before it escapes the wind to 
encounter the interstellar medium. There are at least four regions of the wind that are relevant 
to the relativistic phase of the afterglow (jRamirez-Ruiz et al.ll2005l ): from the inside out, these 
are an expanding WR wind {po/c^ = M^fj^i/AirvwRr'^, where r is the distance from the star), the 
shocked WR wind (po ^ constant), the shocked RSG wind {po ~ constant), and a freely expanding 
RSG wind {po/c^ = M^r/ ^^ttv^sg'''^) ■ Beyond these lie another shocked part of the RSG wind, the 
shocked ISM, and finally the unshocked ISM. 

Density inhomogeneities in such a stratified structure could be produced by several processes. 
First of all, the acceleration region of the WR wind, which extends to a fe w times the stellar r adius, 
is known to be clumpy. Emission- line data indicate ac celerating "blobs" (jMoffat et al.lll988l ) with 
density filling factors / ~ 0.05 — 0.25 (ICrowtheii l2006l . and references therein). These must occur 
where the continuum is optically thin and must be large enough transversely to cover an appreciable 
part (> 10%) of the stellar photos phere. The c l umpin ess of this region is in accord with theory, 
since the line-driving mechanism of lCastor et al.l (Il975l ). which explains the gross properties of WR 



and main-sequence 0-star winds rather satisfactorily, is known to be unstable both radially and 
nonradially on scales larger than the "Sobolev length" /sob ~ '''^^th/^wR ~ 10~^r ([Pessart fc: Owocki 
20051 . and references therein), where v^^i ~ 20 kms"^ is the thermal velocity in the wind. Clump 
dimensions at least as small as ~ 10^*^ cm with density contrasts ~ /^^ ~ 10 are therefore expected. 
These clumps may dissolve beyond the acceleration region (r ^ 10^^ cm) if the optically thin wind 
maintains a uniform temperature comparable to the color temperature of the star; the clumps are 
then at higher density than their surroundings and will expand on their sound-crossing time. After 
crossing t he reverse shock, the wind reaches temperatures ~ 10® K, a regime that is thermally 
unstable ( Fieldlll965l ). Beyond the contact discontinuity, the shocked RSG wind lies at T ~ 10® K, 
which is even more unstable. Thermal instability may give rise to new clumps, whose minimum 
size is controlled by thermal conduction and is therefore very uncertain because the conduction 
rate is probably sensitive to magnetic field. Furthermore, the contact discontinuity is subject to 
Rayleigh - Taylor instability as the shocked RSG wind is accelerated by the less dense WR wind; 
however, iGarcia-Segura &: Francol (|l996l ) found that clump formation is not efficient for the case of 



- 13 - 



pure Rayleigh- Taylor instabilities in thick shells and concluded that a necessary condition for clump 
development to occur is that the shocked s hell be thin enough to allow Vishniac-like instabilities 
(j Vishniad 1 1 983l ) . iRamirez-Ruiz et alj (|2005l ) point out that rapid cooling of the shocked RSG wind 
is indeed likely to lead to a thin shell. Clump development in the shell is controlled by the Kelvin- 
Helmholtz damping of the coupling between the Rayleigh- Taylor and Vishniac instabilities, while 
the lifetime of small dense clumps, once formed, probably again depends upon conduction. 

In short, there are ample reasons to expect a strongly inhomogeneous medium ahead of the 
GRB forward shock, but the sizes of the smallest clumps are uncertain. We make no attempt to 
estimate these sizes but simply state below the properties (lengthscale, density contrast and number 
density) that the clumps would have to have in order to produce post-shock turbulence capable of 
amplifying the magnetic field sufficiently. 

It is of course essential for our proposed mechanism that a seed field exist in the pre-shock 
medium. Magnetic and thermal energy densities are generally comparable in the ISM, but much 
less is known about the magnetic component of the winds of early-type stars. Since such winds are 
believed to be driven by radiative rather than centrifugal forces, as in late-type main-sequence stars, 
and since such stars do not possess surface convection zones, the magnetic field might be very small 
a priori. However, magnetic effe cts have been invoked to explain X-ray emi s sion and anisotrop v 



in the winds of early-type stars (jWolf et al 



199S 



ud-Doula Owocki 



2002 



Schulz et al.ll2003l ). 



Furthermore, since the amplification of a weak (kinematic) magnetic ene rgy density by turbulence 



proceeds exponentially on the timescale of the energy-bearing eddies (jSchekochihin et al 



2002 



and references therein), the number iVeddy of eddy-turnover times required to reach saturation is 
only logarithmically dependent on the strength of the seed field. Therefore, in the estimates of 
the required turnover time made below, we shall simply assume that the magnetic energy fraction 
ee in the pre-shock medium is comparable to that in the ISM, ~ 10~^, so that the number of 
eddy rotations necessary to explain the value of inferred from afterglow observations is of order 

iVeddy ~ 10. 



4. Results 



The formalism described in §11 can be used to predict the eddy-turnover time and vortical 
energy fraction of the turbulent motions produced by an ultra-relativistic shock sweeping-up a 
clumpy medium; if the pre-shock average {i.e., smoothed over clumps) density profile is known, 
a comparison with afterglow models will then let us constrain the lengthscale, overdensity and 
volume filling factor of the circumburst inhomogeneities. Unfortunately, o bservations have not 
yet clearly determined the density profile ahead of the GRB forward shock ([Chevalier &: Lill200d : 
Panaitescu &: Kumaii l200ll . |2002| : I Chevalier et al.l 12004 ) , so that we will consider both free winds 



(energy density /Oq oc r ^) and shocked unif orm winds {po k, constant ) as possible circumburst 



media, keeping in mind that, as described by iRamirez-Ruiz et al.l (|2005l ) and outlined in §111, the 
actual surrounding environment is much more complex. 
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Assuming that the blast wave is adiabatic and effectively spherical and that Eiso = 10^^£'iso,53 erg 
is its isotropic equivalent energy, as derived from the gamma-ray output, we can compute the decel- 
eration radius of the GRB forward shock when about half of the initial energy has been transferred 
to the shocked matter; for an initial Lorentz factor Tq = lO^ro.2, the typical swept-up mass where 
this happens is 

Mdec = ^ ~ 5.6 X 10-6 Si3o,53 Mq . (26) 

i Q C 

For the two circumburst density profiles mentioned above, the deceleration radius of the shock is 
then 



rde 



Mjs^ ^ 3.5 X 10^5 ^i,„_53 ^;^^^8.3 ^wR -5 2 "11 for po oc r'^ 

(ll^^) ' ' « 1-2 X 10^7 (ii;iso,53/nisM,o)'/' roJ/=^ cm for po ^ const. 



where in the first case we have chosen typical WR wind parameters (Mwa = 10 ^Mwr,,-5 M0yr ^ 
and VyjYi = 2000 

''^wR,8.3 kms ^), whereas in the second case a baryon number density comparable 
to the ISM value (tiism = nisM,o cm"'^) has been assumed. This corresponds to a deceleration time 
in the post-shock rest frame 

= V2— ~ / ^ 10^ ^iso,53 ?^WR,8.3 _5 r(^^ S for po OC r'^ ^^^^ 

Tqc I 5.7 X lO'* (£;iso,53/'raisM,o)^''^r^2'^^ s for Po ~ const, 
since the Lorentz factor of the post-shock material in the pre-shock rest frame is ro/\/2- 



4.1. Clump properties 

We consider a pre-shock medium with number density N^fi of identical clumps. Each clump is 
characterized by a gaussian overdensity profile with central density contrast (5max and typical size 
L; if the origin of the axes is in the center of the clump and R, zq and 9 are pre-shock cylindrical 
coordinates, we choose a clump energy density profile 

P0{R, Zo) = Pext[l + <5maxe-(^'+^o)/^'] , (29) 

where pext is the energy density of the inter-clump homogeneous medium. In order to use the 
formalism introduced in §11 — where the density contrast 6 was defined with respect to the mean 
pre-shock energy density po (averaged over clumps), whereas here the inter-clump density pext has 
been used — we should set 

5iR,zo) = ^[1 + 5^axe-(^'+^o)/^'] - 1 ; (30) 
Po 

however, the ratio between /Oext and po can be easily computed for the density profile in eq. ([29]) : 
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provided that 5max-^c,o^^ ^ 1- In this section, we will always assume Pext/po ^ 1; so that the 
formulae in §11 can be used with a pre-shock overdensity profile 

5{R,z,) = 6^,^e-^'''+^oyL\ (32) 

For small density contrasts, we can perform analytical calculations taking into account just 
the leading (first) order term in 5 within eq. (jlSp : the corresponding vorticity will be referred to 
as the "leading-order vorticity" oJloj and the Hankel transform of its projected vorticity will be 
used in eq. ([25]) to analytically compute the leading (second) order in 5max of the vortical energy 
fraction (evort.Lo)- This approximation is certainly well justified for Sma.^ *C 1, but Figure [2] shows 
that the full numerical computation for e^ort is still in reasonable agreement with our analytical 
approximation evort,Lo even for 5max ~ 1- For larger central overdensities a numerical calculation is 
required, and we could fit the numerical results with a fitting function 

= '™^^'^°l + Ci(Lx)- ^^^^ 



(ci and C2 are fitting parameters), so that for small overdensities (5max ^ 1) we recover the result 
of the analytical computation. 



4.2. Eddy-turnover time 

Vorticity embedded in clumps by the passage of the GRB forward shock can be responsible 
for the magnetic fraction inferred from afterglow models only if vortical motions are fast and 
energetic enough to amplify the field up to the observed value before the shock deceleration time, 
when adiabatic expansion would significantly reduce the particle energies available for the afterglow 
emission. 

In the leading-order approximation described above, an estimate of the eddy-turnover time for 
the density contrast in eq. (f32]) is 

teddy = 1 ~ 6.5 X 10^ Li3-^ s (34) 

|a;Lo(^,0)| (^max 

where L = W^^Li^ cm; the reference value chosen for the central overdensity (5max is in agreement 
with observations (see §111) and reasonably satisfies the requirements for the leading-order approx- 
imation (see Figure [2]). If A'eddy = lOA^^eddy,! is the number of eddy rotations necessary to explain 
the observed ee, the requirement A'eddyteddy < ^dec for a wind-like or ISM- like circumburst medium 
gives respectively 



^13 < 0.3 Ei,a,53 VwR,8.3 ^cddy.l PO OC r ^ 

^133;^ < 8.7 (i?iso,53/nisM,o)'/' r-^/=^ ^ for po « const. 



(35) 



Let us emphasize that, for smaller initial shock Lorentz factors Fq, a larger size L and weaker 
overdensity 6raax are enough to satisfy the constraint A'cddytcddy < *dcc- 
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4.3. Vortical energy fraction 

The kinetic energy density invested in the vortical part of the post-shock flow should be 
a significant fraction of its proper energy density, since — as stated in §11 — we assume that the 
amplified magnetic fraction, which is required to match the observational inference ^ 10^^, 
will eventually be comparable to the vortical fraction, ~ fivort- However, ea might be smaller 
than evort if the backreaction of magnetic field on turbulent motions is important well before 
equipartition between the magnetic and vortical energy densities; on the other hand, evort may also 
be a lower limit for if secondary shocks created by the overlapping sound waves of many different 
overpressured clumps significantly contribute to magnetic field amplification. 

The results of a numerical calculation for the vortical fraction of clumps with overdensity profile 
in eq. (j32p are shown in Figure [3] {solid line), but for small central overdensities the leading-order 
analytical computation described above gives a reliable estimate of the magnetic energy fraction 
produced by turbulence in GRB afterglows: 

6b ~ 6™rt,.o - 1.8 X 10-3 f -1 (^) ' ^ (36) 

where F = 10^ r2 and the reference values for the central overdensity (5max and the clump volume 
filling factor NcflL^ have been chosen in order to match the observational constraint ^ lO^^. So, 
clumps with moderate density contrasts {Smax ~ 2, in agreement with the observational evidences 
reviewed in §111) can justify the lower limit of the magnetic energy fraction inferred from afterglow 
models. Higher density contrasts (see Figure [3]) would be necessary to achieve larger magnetic 
fractions; however, it is worth recalling that our model is applicable only under the assumptions 
that the clump overdensity is not larger than ~ F (see §n) and that (5max Nc,oL^ < 1, so that the 
overdensity profile can actually be described by eq. ([32]) . 

Eq. (|36p suggests that smaller overdensities would be enough to satisfy the observational con- 
straints on ea as the forward shock slows down, since oc f~^; the magnetic fraction is then 
expected to change during the afterglow stage, and in §V we will discuss the possible observational 
consequences of its time evolution. 



5. Summary and Discussion 

We have proposed that the post-shock magnetic fields of GRB afterglows may arise from macro- 
scopic MHD turbulence rather than microscopic plasma instabilities. The source of turbulence is 
vorticity produced when the shock encounters density inhomogeneities in the pre-shock medium. 
We presume that the magnetic energy fraction (ee) that results is comparable to the energy frac- 
tion of the turbulence. The ul tra-relativistic Geometrical Shock Dynamics formalism of Paper I 



(| Goodman &: MacFadyenI 120071 ) allows an easy, though approximate, calculation of the vorticity 



produced by a given density inhomogeneity in the limit that the shock Lorentz factor F ^ 1. 
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In this picture, the observational inference that ^ 10~^ constrains both the ampUtude and 
the lengthscale of inhomogeneities. Eqs. (j33p and (j36p roughly relate the total energy fraction in 
vortical motions to the volume filling factor and density contrast of the clumps; this energy fraction 
must be comparable to the inferred post-shock e^- Filling factors and density contrasts of order 
unity are required when the shock is still highly relativistic. Eq. (jSSp . on the other hand, express 
the constraint on the lengthscale and density contrast of individual clumps (independently of their 
volume filling factor) so that the eddies can wind the magnetic field up to the observed value in 
less than the expansion time of the shock. This second constraint favors small lengthscales, so that 
the clumps responsible for field amplification w ould probably be too s mall, at least individually, 



to modulate the afterglow light curve. Actually, iNakar &: GranotI (|2006l ) have recently shown that 



any pre-shock inhomogeneity can hardly be responsible for the rebrightenings observed in some 
afterglow light curves. 

There is a question whether a fiuid treatment of the post-shock flow is justified at all since 
the plasma is collisionless. The same question arises in supernova remnants, to which the standard 
answer is that magnetization of the plasma ensures a short effective mean free path. The present 
case is more extreme because the particles are relativistic, the lengthscales on which fiuid- like be- 
havior is required are smaller, especially in the present work, and the pre-shock field is energetically 
negligible. The relativistic Larmor radius of the post-shock ions based on the compressed ambi- 
ent field is Tl ~ TnipC^ / eB ~ 10^^(3 /iG/Bq) cm, where Bq is the pre-shock field strength. This 
is smaller than the maximum tolerable clump size for field amplification at the beginning of the 
afterglow phase, though only barely so ( §4.2p . Furthermore, whether or not the Weibel instability 
can produce persistent magnetic fields, it should enforce fluid-like behavior by isotropizing particle 
distribution functions whenever counterstreaming plasmas overlap. 

A basic conclusion of this work is that vortical turbulence becomes easier to produce with 
decreasing shock Lorentz factor. Both the energy and timescale constraints become easier to satisfy 
as r decreases (however, our ultra-relativistic approximations break down as F ^ 1). Therefore, if 
the post-shock magnetic energy density is produced by macroscopic turbulence, it is likely that ee 
will evolve as the shock ages, complicating the task of drawing physical inferences about the GRB 
environment from the observational data. The abundance of early X-ray light curves provided by 
the Swift satellite has already led to models that are m ore complicated tha n the rather simple 



theoretical description of the sparser BeppoSAX results (jGalama et al.lll998l ). The light curves 



are not single power laws in time, but show br eaks and sometime s "flares", suggesting a need for 



extended energy input from the central source (jZhang et al.ll2006l . and references therein). But, to 



date, most modelers have assumed constant and eg (the post-shock energy fraction in relativistic 
e^) within individual afterglows, although these parameters are often allowed to vary from one 
afterglow to another. 

The effect of an evolving depends upon the relationship between the observed frequency 
and certain critical frequencies in the synchrotron spectrum. Particularly important is the cooling 
frequency, the doppler-shifted synchrotron frequency of a post-shock electron or positron that 
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radiates much of its energy on a timescale comparable to the age of the shock. For a pre-shock 
medium with mass density profile Pm{f) = c~'^po{r) = Kr~'^ averaged over clumps, and for an 
adiabatic relativistic shock with constant isotropically equivalent energy E^iso (notwithstanding the 
above-cited inferences from Swift data), the cooling frequency evolves as (up to dimensionless 
constants of order unity) 



~ ^"^e \ - ( M-3/2 - 



(1 + 2;grb)o-? 



1 + Zgj^b) t 



(3w-4)/2(4-w) 

(37) 



where rg is the shock radius, t the astronomical observer's time and Zgrb the GRB cosmological 
redshift. Evidently, the shock energy and Lorentz factor scale out of the cooling frequency when the 
latter is expressed in terms of the shock radius. If one could be confident that the early afterglow 
evolves in a freely expanding wind (u; = 2), which seems a priori likely in collapsar models, then 
the evolution of ee could be inferred by measuring that of i^cooi- 

Present evidence suggests that the cooling frequency lies below the X-ray regime in the early 
afterglow phase. This conclusion rests on the usual assumption that the synchrotron-emitting 
electrons are injected with a power law distribution of energies, N{'y)d'y oc ^~^d'y for 7 > ^min ^ l, 
in which p > 2 so that the total energy is dominated by electron energies near 7min'7ieC^i whose 
characteristic observed frequency is i^min- The observed specific flux is often described as a power 
law in time and frequency, F,^ oc t~°'u~^, despite various breaks and the aforementioned flares. 
If synchrotron emission dominates and fcool > i^mm (slow-cooling regime), the spectral index is 
[3 = (p — l)/2 if z^min < v < i^cooi and (3 = p/2 \i v > fcooi- It is believed that the acceleration index 
p is not much larger than 2, perha ps p k, 2.2 — 2.3, w hich is then consistent with the observed X-ray 



indices /3 ~ — 1 observed by Swift (jZhang et al.ll2006l . and references therein) only if /it'cooi < 1 keV 



Evaluating eq. (I37p at the deceleration radius appropriate for a Wolf-Rayet wind (see eq. (I27p ). we 
find that the afterglow phase begins with a cooling frequency that is plausibly consistent with this 
constraint: 

/i^^cool.dec ~ 0.2 (1 + Zc.RB)-Hl0^eB)-'/' ^^w^S.S ^WR,-5 ^-o-53 keV. (38) 

Unless the cooling frequency passes through the observed band, one cannot learn much about the 
evolution of 63 from observations at > Vcooi > i^min, because (up to numerical factors of order 
unity) 

47r4 (1 + Zgrb)-i t V7n 

in this regime, where cIl is the luminosity distance; enters the above expression only via the 
correspondence « (1 + Zgrb)~^ (l^)'^ i^BPmy^'^ e/nT'e between observed frequency u and electron 
Lorentz factor 7, and hence is raised to the small exponent {p — 2)/4 < 0.1. On the other hand, 
relation (j39p indicates that the flux a bove the cooling frequency p rovides an excellent measure of 



^i^^^ . . " w.J r^) (39) 



the energy in the electron population (jFreedman &: Waxman 



2001 



-19- 



We thank Ed Jenkins and Bruce Draine for discussions of small-scale inhomogeneities in the 
ISM, and Eli Waxman for discussions of recent GRB afterglow observations and their implications. 
This work was supported in part by the NSF Center for Magnetic Self-Organization at Princeton. 



-20- 



REFERENCES 

Castor, J. I., Abbott, D. C, k Klein, R. 1. 1975, ApJ, 195, 157 

Chevalier, R. A., Fransson, C, & Nymark, T. K. 2006, ApJ, 641, 1029 

Chevalier, R. A., & Li, Z.-Y. 2000, ApJ, 536, 195 

Chevalier, R. A., Li, Z.-Y., & Fransson, C. 2004, ApJ, 606, 369 

Coburn, W., k Boggs, S. E. 2003, Nature, 423, 415 

Crowther, P. A. 2006, ArXiv Astrophysics e-prints 

Dessart, L., & Owocki, S. P. 2005, A&A, 437, 657 

Eshraghi, H. 2003, Physics of Plasmas, 10, 3577 

Field, G. B. 1965, ApJ, 142, 531 

Frederiksen, J. T., Hededal, C. B., Haugb0lle, T., & Nordlund, A. 2004, ApJ, 608, L13 
Freedman, D. L., & Waxman, E. 2001, Astrophys. J., 547, 922 

Galama, T. J., Wijcrs, R. A. M. J., Bremer, M., Groot, P. J., Strom, R. G., Kouveliotou, C, & 
van Paradijs, J. 1998, ApJ, 500, L97+ 

Garcia-Segura, G., & Franco, J. 1996, ApJ, 469, 171 

Genet, F., Daigne, F., & Mochkovitch, R. 2007, ArXiv Astrophysics e-prints 
Goodman, J., & MacFadyen, A. 2007, J. Fluid. Mech., submitted (Paper I) 
Ishizuka, T., Hashimoto, Y., & Ono, Y. 1964, Progress of Theoretical Physics, 32, 207 
Lazzati, D., Rossi, E., Covino, S., GhiseUini, G., & Malesani, D. 2002, A&A, 396, L5 
Medvedev, M. V., k Loeb, A. 1999, ApJ, 526, 697 
Meszaros, P., & Rees, M. J. 1997, ApJ, 476, 232 

Moffat, A. F. J., Drisscn, L., Lamontagne, R., & Robert, C. 1988, ApJ, 334, 1038 
Nakar, E., & Granot, J. 2006, ArXiv Astrophysics e-prints 
Nakar, E., Piran, T., & Granot, J. 2003, New Astronomy, 8, 495 
Panaitescu, A. 2005, MNRAS, 363, 1409 
Panaitescu, A., & Kumar, P. 2001, ApJ, 554, 667 



- 21 - 



— . 2002, ApJ, 571, 779 

Ramirez-Ruiz, E., Garci'a-Scgura, G., Salmonson, J. D., &; Perez-Rendon, B. 2005, ApJ, 631, 435 
Sari, R. 2006, Physics of Fluids, 18, 027106 

Schaefer, B. E., Gerardy, C. L., Hoflich, P., Panaitescu, A., Quimby, R., Mader, J., Hill, G. J., 
Kumar, P., Whcclcr, J. C., Eraclcous, M., Sigurdsson, S., Meszaros, P., Zhang, B., Wang, 
L., Hessman, F. V., & Petrosian, V. 2003, ApJ, 588, 387 

Schekochihin, A. A., Cowley, S. C., Hammett, G. W., Maron, J. L., & McWilliams, J. C. 2002, 
New Journal of Physics, 4, 84 

Schulz, N. S., Canizares, C., Huenemoerder, D., & Tibbets, K. 2003, ApJ, 595, 365 

Silva, L. O., Fonseca, R. A., Tonge, J. W., Dawson, J. M., Mori, W. B., & Medvedev, M. V. 2003, 
ApJ, 596, L121 

Spitkovsky, A. 2005, in AIP Conf. Proc. 801: Astrophysical Sources of High Energy Particles and 
Radiation, ed. T. Bulik, B. Rudak, & G. Madejski, 345-350 

ud-Doula, A., & Owocki, S. P. 2002, ApJ, 576, 413 

Uhm, Z. L., & Beloborodov, A. M. 2007, ArXiv Astrophysics e-prints 

Vishniac, E. T. 1983, ApJ, 274, 152 

Volk, H. J., Berezhko, E. G., & Kscnofontov, L. T. 2005, A&A, 433, 229 
Wang, X., & Loeb, A. 2000, ApJ, 535, 788 

Waxman, E. 2006, Plasma Physics and Controlled Fusion, 48, B137 
Whitham, G. B. 1974, Linear and Nonlinear Waves (Wiley) 

Wolf, B., Stahl, O., k Fullerton, A. W., eds. 1999, Variable and Non-spherical Stellar Winds in 
Luminous Hot Stars 

Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 

Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 459 

Zhang, B., Fan, Y. Z., Dyks, J., Kobayashi, S., Meszaros, P., Burrows, D. N., Nousek, J. A., & 
Gehrels, N. 2006, ApJ, 642, 354 

Zhang, B., Kobayashi, S., & Meszaros, P. 2003, ApJ, 595, 950 



This preprint was prepared with the AAS lAT^^X macros v5.2. 



- 22 - 




(c) i (d) 

Fig. 1. — History of a density clump overrun by an ultra-relativistic shock (mean Lorcntz factor T), 
viewed in the shock rest frame, (a) The Lorentz-contracted (by a factor ~ f spherical clump 
approaches the shock (thick line) with v ^ c. (b) After the shock, the clump contracts further by 
xl/3 (or x(2\/2)~^ in its rest frame) and moves downstream at c/3. (c) The clump re-expands 
to reach pressure equilibrium and emits sound waves or weak shocks, (d) Vorticity created by the 
shock passage begins to roll up the clump. 
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Fig. 2. — Ratio between the exact numerical calculation of the vortical energy fraction (cvort) and 
its analytical approximation (evort.Lo in eq. (f36]) . keeping only the leading order in the central over- 
density (5max within the vorticity u) as a function of ^max! deviations from the anaytical calculation 
are significant only for (5max ^ 1- Dotted line: Taking into account, in the numerical calculation, 
only the non-integral contribution to the vorticity (i.e., neglecting the second row in eq. (jlSp ). Solid 
line: Full numerical calculation, including both the integral and non- integral terms in eq. ()15|) : the 
non-integral contribution is found to be dominant even for large overdensities. 
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Fig. 3. — Solid line: Full numerical calculation of the vortical energy fraction evort as a function 
of the central overdensity 5max, assuming an average shock Lorentz factor F = 100 and a clump 
volume filling factor NcflL^ = 0.25 (a different choice for this parameter would simply shift the 
curve). Dotted line: Analytical best-fit with the fitting function in eq. and fitting parameters 
ci ~ 0.176 and C2 ~ 1.054. The area between the dashed lines gives the region allowed by afterglow 
observations for the magnetic energy fraction eg. 



